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ABSTRACT 

NASA’s Genesis sample return mission will be the first to return material from beyond the 
Earth-Moon system. NASA Langley Research Center supported this mission with aerothermody- 
namic analyses of the sample return capsule. This paper provides an overview of that effort. The 
capsule is attached through its forebody to the spacecraft bus. When the attachment is severed 
prior to Earth entry, forebody cavities remain. The presence of these cavities could dramatically 
increase the heating environment in their vicinity and downstream. A combination of CFD calcu- 
lations and wind tunnel phosphor thermography tests were employed to address this issue. These 
results quantify the heating environment in and around the cavities, and were a factor in the deci- 
sion to switch forebody heat shield materials. A transition map is developed which predicts that 
the flow aft of the penetrations will still be laminar at the peak heating point of the trajectory. As 
the vehicle continues along the trajectory to the peak dynamic pressure point, fully turbulent flow 
aft of the penetrations could occur. The integrated heat load calculations show that a heat shield 
sized to the stagnation point levels will be adequate for the predicted environment aft of the pene- 
trations. 


NOMENCLATURE 

d depth of cavity [in, cm] 

H altitude [km] 

h heat transfer coefficient [kg/nr-s] 

Q heat load [J/cnr] 

q heating rate [W/cm : ] 

R radius [in, m] 

Re Reynolds number 

Re tv „ cell Reynolds number at wall 

T temperature [K] 

t time [,vj 

V velocity [m/s] 

w’ width of cavity [in, cm] 

x radial distance from spin axis [in, m] 
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z axial distance from nose [in, m] 

a angle of attack [deg] 

P forebody cone half angle [deg] 

8 boundary layer thickness [in. m] 

e surface emissivity 

0 momentum thickness [in, m] 

0 angular location of cavity [deg] 

p density [kg/nr] 

1 cavity taper angle [deg] 

Subscripts 


b base 

c cavity 

D maximum diameter of model 

e boundary layer edge 

FR Fay-Riddell formula 

lam laminar 

n nose 

nom nominal cavity 

s shoulder 

turb turbulent 

w j wall 

<*> free stream 

0 momentum thickness 
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INTRODUCTION 

The Genesis mission is part of NASA's Discovery 
Program, which investigates answers to the fundamen- 
tal questions of the solar system. The objective of 
Genesis (http://genesismission.jplnasa.gov/) is to col- 
lect solar wind from the Sun in order to gather informa- 
tion about the origins of the solar system. These are to 
be the first samples brought back to Earth from beyond 
the Earth-Moon system. 

Set to launch in mid-2001, the Genesis spacecraft 
(shown in Figure 1) will cruise for three months to a 
halo orbit about the LI libration point (0.01 AU from 
Earth). There it will deploy collector arrays which will 
be exposed to solar wind for two years. At the end of 
that time, the collectors will be stowed in a sample re- 
turn capsule (SRC), and the spacecraft will return to the 
Earth. Upon its arrival in 2004, the SRC will separate 
from the rest of the spacecraft and reenter the atmos- 
phere, decelerating with the aid of a parachute. A heli- 
copter will perform a mid-air recovery over the United 
States Air Force (USAF) Utah Test and Training Range 
(UTTR). 

A drawing of the Genesis SRC geometry is pre- 
sented in Figure 2, and its defining parameters are given 
in Table 1. Its forebody (a spherically-blunted cone 
with a rounded shoulder) is similar to that of the Star- 



Figure 1. Genesis spacecraft with collector arrays 
deployed. 



Figure 2. SRC layout. 


Table 1. Parameters for SRC Geometry. 


Parameter 

Baseline 

Updated 

P [deg] 

59.81 

59.81 

R„[m] 

0.4394 

0.4298 

R h [m] 

0.7448 

0.7543 

R< t«i] 

0.0332 

0.0335 


dust 1 5 SRC. As with Stardust, NASA Langley Research 
Center (LaRC) has provided support to Lockheed- 
Maitin Astronautics in their design of the Genesis SRC 
through entry aerodynamic, aerothermodynamic. and 
flight dynamics 6 ' 8 analyses of the capsule. This paper 
provides an overview of the aerothermodynamics ef- 
fort. 


The design of the Genesis spacecraft requires that 
the SRC be attached to the rest of the vehicle via three 
structural inserts of non-TPS material which pass 
through holes in the forebody heat shield. In this way, 
the collectors can be deployed from the aft side of the 
SRC. Figure 3 shows the size and locations of these 
heat shield penetrations. When the SRC separates from 
the spacecraft bus, the bolts are severed and retracted to 
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Cavity Cross Section 


Figure 3. Penetration geometry on the forebody of 
the SRC flight article. 
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form the cavities. The presence of these cavities during 
reentry may result in high, localized heating at the 
downstream edge of the cavities due to flow separation 
and reattachment within the cavities, accompanied by a 
heating augmentation downstream of the cavities. 
These penetrations could also induce transition to tur- 
bulence. which in turn could dramatically increase the 
heating environment aft of the penetrations. 

Wind tunnels tests for the Apollo Command Mod- 
ule 9 and the Viking 10 aeroshell show large increases in 
heating (greater than three times stagnation levels) im- 
mediately downstream of penetrations. Further down- 
stream, heating levels were observed to be consistent 
with the onset of transition to turbulent flow r . It could 
also be seen that the downstream heating levels are a 
function of cavity size. 

The above results are not directly applicable to the 
Genesis SRC (due to differences in capsule geometry, 
cavity geometry, etc.), so computational fluid dynamics 
(CFD) calculations and wind tunnel phosphor thermo- 
graphy tests were employed. Since the cavities are very 
small at the scale of a wind tunnel model, quantification 
of the magnitude of the heat pulse within, and in the 
immediate vicinity of, the cavities was handled com- 
putationally. Wind tunnel tests using the phosphor 
thermography technique were used to determine the 
heating augmentation aft of the penetrations for a ma- 
trix of freestream conditions and cavity geometries. 
Finally, these experimental results were extrapolated to 
tlight conditions using CFD. These results were a factor 
in the decision by Lockheed- Martin Astronautics to 
switch from the baseline PICA (Phenolic Impregnated 
Carbon Ablator) heat shield material to carbon-carbon. 

ANALYSIS 

Numerical Approach 

Laminar, non-blowing, thin-layer Navier-Stokes 
solutions are generated with the Langley Aerothermo- 
dynamic Upwind Relaxation Algorithm (LAURA). 11 
LAURA is a finite volume, shock capturing algorithm 
capable of calculating flows in chemical and thermal 
nonequilibrium. It can be used for inviscid, thin-layer 
Navier-Stokes, or full Navier-Stokes computations. The 
solutions presented herein are for the thin-layer Navier- 
Stokes equations, except for one full Navier-Stokes 
solution in the sensitivity study. Turbulence calcula- 
tions are performed using the Baldwin-Lomax 12 alge- 
braic model. 


For the flight cases, the flow was assumed to be in 
chemical (11-species air model) and thermal (separate 
translational and vibrational temperatures) nonequili- 
brium. For the surface catalysis, mass fractions were set 
to their values at the freestream temperature, resulting 
in a conservative estimate of the heating rate. Binary 
diffusion fluxes were corrected so that their sum equals 
zero. For the Bight cases, one of two approaches was 
used: I) T w = 2900 K is specified, or 2) surface wall 
temperatures were computed using the radiative equi- 
librium wall temperature assumption with an assumed 
emissivity of 0.9. 

For the wind tunnel cases, the flow was treated as a 
perfect gas. A constant wall temperature of T w = 300 K 
was specified. 

The numerical results presented herein were com- 
puted on SGI machines (32-bit word length). A single 
case was repeated on a Cray machine (64-bit word 
length) to verify that results are consistent. In all cases, 
solutions were judged to be converged w hen the surface 
heating values at selected points over the vehicle 
changed less than 0.01 per cent over 1000 iterations. 
Typically, 20.000 iterations were required to reach con- 
vergence. 

Experimental Method 
20-lnch Mach 6 Tunnel 

With its high freestream Reynolds number range 
(1 .6 to 29.3 x 10 6 m y ), the LaRC 20- Inch Mach 6 facil- 
ity is well suited for looking at transitional and turbu- 
lent flows. The tunnel has a two-dimensional contoured 
nozzle, with the bottom and top walls being contoured 
and the sidewalls parallel. The two-dimensional nozzle 
provides very uniform flow at the test section. Tem- 
pered glass windows on the top and two sides allow for 
a significant amount of visual access for the phosphor 
systems. Models are mounted on an injection system in 
a housing below- the test section. During a run, the 
models can be injected in as little as 0.5 s (out of a total 
tunnel run time of 4 s). 

Phosphor Thermography Technique 

The two-color phosphor thermography technique, 
as developed at NASA LaRC. n was used to obtain sur- 
face heat transfer data in the wind tunnel. With this 
method, ceramic wind tunnel models are coated with 
phosphor crystals, which fluoresce (when illuminated 
with ultraviolet light) with an intensity that is tempera- 
ture dependent. 
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During a wind tunnel run, intensity images are ob- 
tained with a state-of-the-art image acquisition system. 
Afterwards, the resulting intensity images are converted 
to quantitative temperature mappings using a weighted 
two-color relative-intensity fluorescence theory. Tem- 
perature data can be acquired even on wind tunnel 
models with complex geometries. 

Once quantitative temperature mappings have been 
obtained, the next step is to derive heat transfer map- 
pings. Given surface temperature variations over time at 
a location on a model, the heat transfer rate at that loca- 
tion can be inferred from one-dimensional heat con- 
duction equations. Global heat transfer data can be effi- 
ciently obtained within minutes of a wind tunnel run. 
The LaRC phosphor thermography methods allows a 
three-orders-of-magnitude increase in the amount of 
heat transfer data gleaned from a given wind tunnel run, 
with an order-of-magnitude reduction in cost and time, 
as compared to conventional discrete gauge methods. 
This allows aeroheating information to be brought for- 
ward earlier in the design process than was previously 
possible. 

The initial temperature and run temperature are the 
dominant uncertainties when lower surface tempera- 
tures are measured. Indeed, as the surface temperature 
approaches the initial temperature, the uncertainty 
climbs to infinity. The high precision errors in tem- 
perature are responsible for the large amount of scatter 
observed in phosphor heating data obtained from lower 
surface temperature readings. In most cases, though, the 
wall temperature measurements are higher, and in the 
case of the bias and total uncertainties, the uncertainty 
due to thermal properties is the most significant error 
source. For the majority of the temperature range, the 
resulting bias, precision, and total uncertainties are 7-8 
per cent, 3-8 percent, and 8-10 per cent, respectively. 

Model Description 

Four different scale (5, 6.67, 8.33 and 10 per cent) 
configurations were fabricated. The resultant model 
maximum diameters are 3, 4, 5, and 6 in (0.0752, 
0.1006, 0.1392, and 0.1504 m), respectively. Four pat- 
terns, one of each desired model scale, were created 
using the stereo- lithography process. Using these pat- 
terns, models were cast with a fused silica ceramic in- 
vestment casting method and coated with phosphor 
crystals. Eight models were fabricated, two of each 
model scale. The designations are as follows: the 3- in 
models are 3A and 3B, the 4-/>r models are 4A/4B. the 
5 -in models are 5A/5B. and the 6 -in models are 6A/6B. 


Each model was manufactured with six cavities (of 
different sizes, radial positions, and depths) on the 
forebody. The cavity specifications are listed in Table 2 
and shown in Figure 4. Cavity angular locations (0) 
start from the left of the figure and proceed in a 
clockwise manner at 60 -deg intervals. The cavities were 
placed at two different radial locations: 40 and 70 per 
cent of the maximum radius (R b ) of the configuration. 
Three different cavity widths (nominally 0.5, 1.0. and 
2.0 in at full scale) are included at the appropriate 
scaling. Nominal cavity depths are 0.03 in, with the 
exception of the cavity at 9 = 30 Ueg which has a 
depth of 0.06 in. 

Table 2. Cavity specifications for 1/10-scale model. 


Angular 
location 
(9) [<%] 

Radial 

location 

(x/Rh) 

w [in] 

d [in] 

0 

0.70 

0.10 

0.03 

60 

0.40 

0.20 

0.03 

120 

0.40 

0.10 

0.03 

180 

0.70 

0.05 

0.03 

240 

0.70 

0.20 

0.03 

300 

0.70 

0.10 

0.06 



Figure 4, Schematic diagram of the cavity 
configuration for the wind tunnel models. 

The cavity axis is normal to the surface. The wall 
of each cavity is canted at x =35 deg away from that 
axis (see Figure 2). In some cases (cavity at 0 = 300 deg 
on model 4B, cavity at 0 = 0 deg on 6B), the cavity (as 
cast) was not satisfactory on the final model. The cavity 
was restored at these locations by drilling a hole normal 
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to the surface to the desired depth. This approach does 
not preserve the 3 5 -deg taper for these cavities, however. 

RESULTS AND DISCUSSION 
Computational Trajectory Focal Points 

Initially, solutions over the forebody were calcu- 
lated for three points (t = 260, 270, 280 s) along the 
nominal trajectory. The stagnation point heat pulse 
(evaluated using the Fay-Riddell relation) and 
freestream Re history from that trajectory are shown in 
Figure 5, along with all trajectory points at which 
NASA LaRC performed calculations. Freestream con- 
ditions for those points are given in Table 3. 



t [S] 


Figure 5. Reference heating and Reynolds number 
history from nominal trajectory. 

Computational Sensitivity Study 

The effect of grid resolution and thin-layer Navier- 
Stokes assumption on the laminar, non-blowing, surface 
heating for the smooth forebody of the Genesis Sample 
Return Capsule (SRC) was examined. LAURA solu- 
tions were generated at the peak heating point of a 
nominal trajectory with grids of varying resolution in 


Table 3. Freestream Conditions at Calculation 
Points Along Nominal Trajectory 


t M 

V (w/v| 

P \kg/m'\ 

T \K\ 

Re | Kl'm'l 

Note 

260 

1049(1 

1. 0230-4 

216.8 

0.0755 


266 

9955 

2.353e-4 

234 4 

0 1080 


270 

9360 

3.78 le-4 

2446 

0.2254 

Peak heating 

276 

8098 

6.920c-4 

255.7 

0.2410 


280 

7092 

9. 7 28c -4 

262.1 

0.4156 

Peak pressure 

284 

6059 

l.312e-3 

264.6 

0.3325 


290 

4625 

1.93 2c -3 

261.4 

0.3775 



both the streamwise and normal directions. Freestream 
conditions at the peak heating point along the trajectory 
(see Table 3) were used for all computational solutions. 
To simplify the calculations, a constant wall tempera- 
ture of 2900 K was assumed. This value falls within the 
range of calculated radiative equilibrium surface tem- 
peratures. 

A grid of 100 streamwise by 64 body-normal cells 
over the baseline SRC forebody was generated for pre 
liminary aerothermodynamic computations. The grid 
contains 30 cells in the streamwise direction on the 
spherical nose cap, 40 cells on the conical flank, and 30 
cells on the rounded shoulder. This grid served as the 
baseline grid for all computations over the smooth fore- 
body. 

In order to assess the effects of grid density and 
clustering on the forebody surface heating, a number of 
laminar solutions were computed over grids derived 
from the baseline grid. The streamwise grid density was 
varied by coarsening the baseline grid in each region 
independently. The grid clustering in the normal direc- 
tion was altered by varying the cell Reynolds number 
(^ e cdl) based on the local speed of sound for the first 
cell away from the wall. The number of cells in the 
normal direction was also varied. Table 4 lists the ma- 
trix of grids. 


Table 4. Grid Resolution Run Matrix 


Grid 

Streamwise Cells 

Normal Direction 

Note 

Nose 

Cone 

Shoulder 

Cells 

Re cdl 

0 

15 

20 

15 

64 

10 

Current LAURA default 

1 

15 

20 

15 

64 

2 

Nettelhorst and Mitcheltree 14 

2 

15 

20 

15 

128 

1 

Fine grid (normal) 

3 

6 

20 

10 

64 

10 

Extra-coarse nose, coarse shoulder 

4 

30 

40 

30 

64 

10 

Fine grid (streamwise) 

5 

10 

20 

10 

64 

10 

Coarse nose, coarse shoulder 
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Effect of Normal Grid Clustering 

Three grids were used to estimate the effect of grid 
clustering in the normal direction on surface heating. 
Referring to Table 4, these are the low wall-resolution 
default grid with Re wU = 10 (Grid 0), a medium wall- 
resolution grid recommended by Nettelhorst and 
Mitcheltree 14 (Grid 1) with Re LvM = 2, and a fine resolu- 
tion grid with 128 cells in the body-normal direction 
and Re c ,n = KGrid 2). The surface heating computed 
over these grids is shown in Figure 6. 



x[m) 

Figure 6. Effect of normal grid clustering on surface 
heating. 

Results from all three grids show similar heating 
levels in the stagnation region with variations at the 
stagnation point of 1 per cent. The primary effect of the 
grid clustering occurs on the conical flank and on the 
shoulder. Using the fine resolution grid (Grid 2) as the 
baseline, the medium resolution grid (Grid 1 ) predicts 
heating levels that are 9 per cent higher at the end of the 
cone and on the rounded shoulder. The low-resolution 
grid (Grid 0) does a better job by predicting heating 
levels that are within 4 and 1 per cent of the baseline 
levels on the cone and shoulder, respectively. The ob- 
servation that the grid with a lower resolution at the 
wall (Re^n = 10) gives heating rates closer to highest 
resolution (Re tl . n = 1) than the grid with the medium 
resolution (Re C(?1 , = 2) appears to be counterintuitive. 
However, for a given number of normal cells, reducing 
the resolution at the wall decreases the amount of grid 
stretching and gives a better distribution of grid points 
across the entire boundary layer. Based on this com- 
parison, 64 cells in the normal direction with Re u , t | = 10 
adequately resolves the laminar heating and will be 
used for the remaining solutions. 


Effect of Streamwise Grid Density 

Four grids with varying numbers of streamwise 
cells were used to estimate the effect of streamwise grid 
density on surface heating. In the nose region, solutions 
were computed using grids with 30, 15, 10. and 6 cells 
(corresponding to Grids 4, 0, 5, and 3 in Table 4, re- 
spectively). Figure 7 depicts the heating levels in the 
nose region for each grid. The heating levels using all 
four grids are very similar with only the 6-cell run 
showing a relatively small deviation of 4 per cent at the 
stagnation point. Since the 10 nose-cell solution is on 
par with the higher resolution results, and the reduction 
in total computational time is minimal for the 6 nose- 
cell grid over the 10 nose-cell grid, 10 cells in the nose 
region will be used for the remainder of this study. 



x [m] 

Figure 7. Effect of streamwise grid density on nose 
heating. 

The effect of the number of streamwise cells on the 
heating levels in the shoulder region is shown in Fig- 
ure 8. Solutions are computed using grids with 30, 15, 
and 10 cells in the shoulder region (corresponding to 
Grids 4, 0, and 5 in Table 4, respectively). The peak 
heating rate on the shoulder predicted by the three grids 
agree to within 3 per cent and occurs at the end of the 
conical flank and beginning of the shoulder region. 1 his 
comparison suggests 10 cells in the shoulder region are 
sufficient to resolve the area around the cone-shoulder 
interface where the local peak in heating occurs, and 
that resolution will be used for this study. 

Results for Grids 4 and 0 (with 40 and 20 cells on 
the cone, respectively) have heating rates along the 
cone within 1 per cent of each other (not shown). Al- 
though an even coarser grid on the flank might ade- 
quately predict heating levels, 20 cells will be used in 
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anticipation of including cavities on the forebody in 
future computations. 



Figure 8. Effect of streamwise grid density on 
shoulder heating. 

Effect of Navier-Stokes Equation Set 

Although not shown, the thin-layer Navier-Stokes 
solution on Grid 0 was repeated using the complete 
Navier-Stokes equation set. There is no discernible dif- 
ference between the two solutions, suggesting that the 
use of the thin-layer Navier-Stokes equations is appro- 
priate for computations over the fore body. 

Laminar Heating Calculations 

Smooth Forebody (No Cavities) 

Laminar surface heating rates over the smooth 
forebody of the baseline SRC were calculated for three 
points (t = 260. 270, 280 ,y) along a nominal trajectory. 
With the sensitivity study as a guide, each solution used 
a grid with 40 cells in the streamwise direction (10 on 
the nose. 20 on the conical flank, and 10 on the rounded 
shoulder) and 64 cells normal to the body. The grid was 
clustered in the normal direction using a cell Reynolds 
number of 10. A radiative equilibrium wall temperature 
w r as prescribed. The surface heating and temperature 
distributions for these cases, which exhibit the expected 
levels as the vehicle traverses the heat pulse, are shown 
in Figure 9. The distribution at peak heating will serve 
as a reference for later calculations that include the ef- 
fects of the cavity, as well as turbulence. 

Axisymmetric Solutions With Cavity 

The impact of the cavity was first addressed by 
modifying the baseline grid to include a dimple. Note 



x[m] 

Figure 9. Surface heating and temperature 
distributions for points along nominal trajectory. 

that this actually models the 3-D cavity as an axisym- 
metric "groove". A cavity with dimensions 0.75 x 0.3 
in (w = 1 .905 cm, d = 0.762 cm) and taper of x = 45 cleg 
was modeled. The conical portion of the baseline grid 
from the earlier sensitivity study was modified to in- 
clude the feature. The resultant grid again has 10 cells 
on the nose and 10 cells on the shoulder. There are 13 
cells between the nose and dimple, 12 cells within the 
dimple (uniformly-spaced), and 16 cells between the 
dimple and shoulder, yielding a grid of dimensions 61 x 
64 cells. 

Axisymmetric calculations were performed for the 
freestream conditions at peak heating (see Table 3 for 
t= 270 v). A constant wall temperature of 2900 K. 
which falls within the range of radiative equilibrium 
surface temperatures seen in Figure 9, was prescribed 
Results in Figure 10 show a spike on the downstream 
lip of the cavity whose magnitude is more than three 
times the stagnation value of 510 W/cm : (see Figure 9). 
The smooth forebody solution is included for reference. 
The surface with the cavity geometry is shown in gray 
(the freestream velocity vector points down and the 
downstream lip is on the right). 

A grid sensitivity study was conducted for the cav- 
ity, and those results are also presented in Figure 10. 
First, the grid density within the dimple was doubled to 
24 cells, and the resolution between the cavity and the 
shoulder increased to 31 cells, yielding a grid with di- 
mensions 88 x 64. A third case increased the dimple 
refinement to 48 cells, and the downstream refinement 
to 43 cells, yielding a 124 x 64 grid. The results show 
that the two highest grid densities give similar results 
for the How separation at the upstream lip However, 
both the breadth and magnitude of the pulse is strongly 


7 

American Institute of Aeronautics and Astronautics 





x [m] 


Figure 10. Effect of grid resolution and geometry on 
heating in vicinity of cavity. 

dependent on the resolution provided by the grid. Still, 
the heating levels of the spike for the 88- and 124-cell 
grids differ by less than 4 per cent. 

With the maturation of the SRC design, a switch 
was made to a cavity with dimensions 1.0 x 0.3 in 
(w = 2.540 cm. d = 0.762 cm) and taper of x = 35 deg. 
The 88 x 64 grid from Figure 10 was modified to match 
the new geometry while maintaining the streamwise 
spacing near the lip. Outside the cavity, the surface grid 
is virtually identical to the 88 x 64 grid. The result is a 
grid with dimensions 89 x 64. Figure 1 1 shows the im- 
pact of that switch on the heating augmentation. Note 
that the primary change (an increase in the magnitude 
of the pulse relative to the 88 x 64 grid) appears to be 
due to taper angle of the cavity wall relative to the fore- 
body inclination. 

A grid independent solution of heating around a 
sharp comer is virtually unattainable, and in reality any 
comer will have a finite radius. Therefore, the effect of 
corner rounding is examined next. First, a radius of 
R c = 0.33 cm was added to the downstream lip of the 
cavity for the 89 x 64 grid. The new grid has 34 cells 
within the dimple and 34 cells between the dimple and 
the shoulder, yielding a grid with 100 streamwise cells. 
A comparison of those two solutions in Figure 1 1 
shows that rounding the cavity corner significantly re- 
duces the heating spike. Still, even with rounding, the 
strength of the spike on the downstream lip of the cav- 
ity is on the order of 2-3 times the stagnation value. The 
grid sensitivity is shown by coarsening the mesh around 
this corner by a factor of two to produce a 94 x 64 grid. 
The magnitude of the resultant heat spike is 6 per cent 
less than that of the finer grid. 


Figure 11. Effect of edge radius and grid resolution 
on heating in vicinity of cavity. 

The next step is to add a rounded comer to the up- 
stream lip. A new grid was generated over the updated 
geometry (see Table 1) using 10 cells on the nose and 
10 cells on the shoulder. It has 16 cells before the dim- 
ple. 46 cells within the cavity (6 on the upstream radius, 
8 on the upstream edge, 12 along the bottom, 8 on the 
downstream edge, and 12 on the downstream radius), 
and 24 cells downstream of it, yielding a grid of dimen- 
sions 106 x 64 cells. A coarsened version of this grid 
omits every other node within the cavity (dropping it to 
23 cells) to produce a grid of 83 x 64 cells. 

Solutions were calculated over these grids at the 
t = 270 5 flight conditions. A radiative equilibrium wall 
temperature boundary condition was used. The results 
are shown in Figure 12. Note that the rounding of the 
upstream lip has less effect on the heating levels than 



Figure 12. Effect of streamwise grid resolution on 
heating in vicinity of rounded cavity. 
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rounding the downstream lip. As one might expect, the 
coarsening of the grid tends to smear the features within 
the cavity, while the peak value changes by 10 percent. 

3-D Solutions With Cavity 

The above axisymmetric solutions were used in the 
preliminary analysis of the candidate cavity geometries. 
A limited number of solutions were computed for a true 
three-dimensional cavity of dimensions 1.0 x 0.3 in 
(w = 2.540 cm, d = 0.762 cm), with R c = 0.33 cm and 
taper of t = 35 c/eg. In preparation for these solutions, a 
singularity-free grid was generated for the updated 
forebody shape (see Table 1 ) with a cavity staggering 
the symmetry plane. (The singularity-free topology has 
no pole boundaries and thus avoids the associated com- 
putational convergence problems.) 

A solution at wind tunnel conditions was computed 
over the singularity-free grid of dimensions 146x 64 
x 64. For comparison, an axisymmetric solution using 
the 100 x 64 cell grid from above was computed at the 
same conditions. Although not shown, the symmetry 
plane solution from the singularity-free grid compares 
well with the axisymmetric solution at the nose, along 
the cone, and at the shoulders. As shown in Figure 13, 
heating in the symmetry plane for the three-dimensional 
penetration is 10 per cent less than axisymmetric dim- 
ple, apparently due to a relieving effect. Note that the 
disconnects in the 3-D solution correspond to the 
boundaries between the 16 computational blocks. The 
breaks are the result of post-processing (when the cell- 


2 5 r 



Figure 13. Laminar centerline heating computations 
in vicinity of cavity at wind tunnel conditions. 


centered results are extrapolated to the surface to pro- 
duce face-centered values, which are then shifted to the 
corner nodes for plotting) and are not indicative of so- 
lution discontinuities. 

Streamline patterns within the cavity are shown in 
Figure 14 for the wind tunnel case. The flow separates 
as it expands around the upstream lip and reattaches on 
the downstream wall. The streamlines are colored by 
local enthalpy, showing the high energy levels after 
flow' reattachment. Contours of local heating for this 
case are shown on the left half of Figure 15 and reveal 
an important piece of information. For the three 
dimensional cavity, the heat pulse maximum is pushed 
off the centerline. That maximum is of the same order 
as the maximum value from the corresponding axi- 
symmetric solution. Thus, contrary to the indications of 
Figure 13, the three-dimensional nature of the flow r 
merely serves to shift the location of the heat pulse 
without significantly changing its intensity. 



Figure 14. Front view of streamlines inside cavity at 
wind tunnel conditions. 

The right side of Figure 15 shows the heating con- 
tours for the peak heating trajectory point. Note that the 
hot spot for flight conditions (maximum q/q, R = 2.95) is 
much stronger than wind tunnel conditions (maximum 
q/q IR = 2.33). and it is located closer to symmetry 
plane. As a result of this intensity at flight conditions, 
the affect of the cavity is felt for a longer distance 
downstream than in the wind tunnel (as indicated by the 
contours at the top of Figure 15). 
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Figure 15. Front view of heating inside cavity at 
wind tunnel and flight conditions. 

Wind Tunnel Results and Repeatability 

The primary objectives of the wind tunnel tests 
were to determine the effectiveness of a cavity as a 
boundary layer trip, to assess the dependence of this 
effectiveness on cavity geometry, to observe the region 
of downstream influence, and determine the level of 
downstream heating augmentation resulting from the 
interaction. A total of 29 wind tunnel runs were con- 
ducted at various Re and orientations. A sampling of 
the resultant phosphor thermography images is pro- 
vided in Figure 16. 

All three images in Figure 16 show augmented 
heating streaks due to the presence of the cavities. Fig- 
ure 16a represents the highest Re tests (Re D = 
3.81 x 1 0 6 ). Heating levels 60 per cent higher than 
stagnation values are observed aft of the cavities. Fig- 
ure 16b is a close-up shot which shows that the heating 
on the downstream edge of the cavity spikes up (as seen 
in the computational results). In figure 16c, the results 
are at a = 15 deg. In this case, the heating aft of the 
cavity has a cusped shape as the flow follows the 
streamlines. (The Genesis SRC will nominally fly at 
a = 0 deg , so the bulk of the results herein are for that 
orientation. Results at angle-of-attack are included in 
the transition map to be presented later. ) 

An assessment of experimental data integrity is 
presented in Figure 17. Figure 17a provides the tempo- 
ral collapse of a sample run at Re D = 1.10 x 10 6 , and 


(a) Model 6A, a=0 deg, Re 0 =3.8I x l(f. 



(b) Model 6A, a—0 deg , Re^—l. 1 1 x 10 6 (I - in cavity). 



(c) Model 6A, a=I5 deg , Re [} = 1.56 x l (f. 

Figure 16. Sample phosphor thermography heat 
transfer images. 
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shows that the heat transfer coefficient levels achieve 
steady-state values. Figure 17b compares runs where 
the only variable is model size, with relatively good 
scalability. Although not shown, a repeat run with the 
same model as well as a run with an identical model 
show’ good (within 3 per cent variability) run-to-run and 
model-to-model repeatability, respectively. 



i.J L X J L L l L I 1. L i 1 L L i-1 1 I L L.1 1 1 1 1 L LL 1 L .i_Ll J.X J_j L. t_L iJ 

0 0 1 02 0 3 04 0 5 0.6 0 7 0 8 0 9 1 0 

x/R b 


(a) Temporal collapse of data for sample run. 



Figure 17. Data integrity test (a = 0 deg). 


Observations from Experiment 

The cavities have two primary effects on the heat- 
ing. First, as seen in the computational results, a large 
heat pulse is typically generated on the down-stream 
wall of the cavity. As can be seen in Figure 16b. the 
presence of this pulse can impact local heating down- 
stream of the cavity. Second, if the flow is sufficiently 
energetic, the presence of the cavity can serve to trip 
the flow' to turbulence (see Figures 1 6a and 16c). 

Figures 18 and 19 present line cuts of the thermo- 
phosphor data from the a = 0 deg runs. In Figure 18. 
each subfigure shows the heating for the six cavities at 
a given Re. The first observation is that the smallest 
cavity produces no discernable heat spike or tripping 
effect at any of the Re tested Further, the How returns 
to laminar levels downstream of all cavities at the low- 
est Re (Figure 18a). For the middle Re (Figure 18b). the 
downstream levels are consistent with a transitioning to 
turbulent flow. At the highest Re (Figure 18c), the flow 
appears to approach fully turbulent values downstream 
of all cavities (save the 0.5-in hole). Thus, for a given 
Re, all but the smallest cavity are equally effective as 
tripping devices. A forward cavity position allows more 
running length for tripped flow to grow to fully turbu- 
lent values downstream. Flowever. for the relatively 
short running length of this forebody, the location of the 
penetration appears to have only a minor impact on the 
heating levels realized aft of the cavity. It should be 
noted that the variation toward the end of the body is 
probably due to the presence of the shoulder (and the 
associated flow expansion). 

In Figure 19, each subfigure shows the heating 
variation with respect to Re for a given cavity. A sub- 
figure for the smallest cavity is not included since its 
influence appeal's to be benign. Note that the subfigures 
for the nominal and 2 -in cavities contain data from both 
radial locations. It is clear from these figures that a 
larger cavity produces the strongest heat pulse. For both 
the nominal and 2 -in cavities, the location nearest the 
stagnation point experiences a stronger heat pulse than 
its downstream counterpart. A comparison of subfig- 
ures 19a and 19b shows little difference between the 1- 
in cavities of depth 0.03 in and 0.06 in. respectively. 
Any differences could be due to the fidelity of the mod- 
eling of the cavities, rather than a dependence on cavity 
depth 
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(a) Re l} = O.SOxICf 




Figure 18. Comparison of phosphor line cuts 
(a = 0 deg), sorted by Rer,. 



(a) Cavity : size: w = I in, d = 0.3 in (0 = 0, 120 deg) 



(b) Cavity size: vv = I in, cl = 0.6 in ( 6 = 300 deg) 



(c) Cavity size: w = 2 in, d - 0.3 in ( 6 =60, 240 deg) 


Figure 19. Comparison of phosphor line cuts 
(a = 0 deg), sorted by cavity size. 
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Figure 20 provides a comparison between experi- 
mental results at Re,, = 1.1 x lit (where the flow ap 
pears to remain laminar downstream of the 0.5-w cav- 
ity) and LAURA computations over a smooth forebody. 
Note the agreement is excellent overall. The largest 
differences are in the region of the shoulder, where the 
phosphor image resolution is insufficient to characterize 
the heating pulse associated with the small radius of 
curvature. 



x/R b 

Figure 20* Comparison of phosphor line cuts with 
computational results (a = 0 deg). 

Transition Map 

A transition map was developed as a means of ex- 
trapolating these wind tunnel data to flight. Computa- 
tional results for the wind tunnel conditions were used 
to calculate boundary layer properties (0, 8) for the 
wind tunnel conditions. These* in turn* were used to 
construct a plot of Re 0 versus w/8 for the experimental 
data set. Figure 21 shows the final product. The set of 
phosphor images were used to define which cavities 
and conditions produced incipient, transitional, or tur- 
bulent flow' downstream. The symbols are colored 
(laminar, turbulent, etc.) according to those interpreta- 
tions of the images. 

The data was augmented with two power law 
curves to estimate the laminar and turbulent boundaries. 
Then a sketch of the flight trajectory (including the 
peak heating and peak dynamic pressure points) was 
overlayed on the figure. Based on the resultant map. the 
flow aft of the penetrations will still be laminar at the 
peak heating point of the nominal flight trajectory. As 
the vehicle continues along the trajectory toward the 
peak dynamic pressure point, a transition to turbulent 


flow aft of the penetrations is predicted Fully turbulent 
tlow could be realized at the peak dynamic pressure 
point. 



Figure 21. Transition map for forebody cavities. 

Note that the numerical solutions show that, for 
wind tunnel conditions, M is subsonic at all of the 
cavities, except on the leeside for runs at angle of at- 
tack. Flight conditions yield supersonic M, at the cavity 
locations. As a result, using Re H rather than Re B /M 0 to 
extrapolate results to flight should give a more conser- 
vative estimate for the heating environment. Additional 
areas of conservatism stem from the fact that wind tun 
nel noise typically trips a flow earlier than in free fight. 

Turbulent Heating Calculations 

The turbulent values were calculated using a 
Baldwin-Lomax algebraic turbulence model. 12 For the 
smooth forebody calculations, a beginning point for 
transition to turbulence was specified at a radial loca- 
tion of x = 0.7 R h . The turbulent solutions were com- 
puted using 80 cells in the normal direction with Re tv)1 = 
2. This gives a better near-wall resolution in the turbu- 
lent region and increases the heating by approximately 
5 per cent over the baseline grid. Note that for the cav- 
ity calculations, the solutions were tripped at the lami- 
nar flow reattachment point on the downstream lip. 

Figure 22 shows laminar and turbulent results at 
wind tunnel conditions (with and without a cavity) 
normalized to the Fay-Riddell stagnation value. For a 
close-up of these laminar results in the vicinity of the 
cavity, see Figure 13. At these conditions, there is in- 
sufficient energy to drive the algebraic model to pro- 
duce fully turbulent flow' when tripped on the smooth 
body. In addition, the turbulent heating levels are quite 
different with and without the cavity at these conditions. 
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Figure 22.Computed laminar and turbulent center- 
line heating distributions at wind tunnel conditions. 

Figure 23 provides a comparison between turbulent 
computations at wind tunnel conditions and the experi- 
mental results. At this Re. the cavity serves as an effec- 
tive trip and the flow downstream of the cavity appeals 
to approach fully turbulent levels. The agreement be- 
tween the CFD results with a cavity and the wind tunnel 
data is not very good downstream of the cavity. This 
would appear to be a problem with the algebraic turbu- 
lence model in the presence of this complicated flow 
feature. Specifically, the ability to resolve the vortical 
flow emanating from the cavity, the applicability of the 
turbulence model (a modified flat-plate relation), and 
the appropriateness of the simple transition model are 
questionable. In reality, the turbulence levels over the 
smooth forebody are more in line with the wind tunnel 
data (see Figure 22 at a similar condition). As a result, 



smooth forebody calculations will be used to estimate 
the turbulent heating at flight conditions. 

Figure 24 presents laminar and turbulent results at 
flight conditions, with and without a cavity. The lami- 
nar results are from Figure 9. At these conditions, the 
flow is sufficiently energized to yield fully turbulent 
flow downstream of the specified trip. Even with the 
cavity, there is sufficient energy to drive the turbulence 
model and produce nearly fully-turbulent heating val- 
ues. Still, the smooth forebody computations provide a 
more conservative value, and will be used in later heat 
load calculations. 


• Laminar, 40 x 64 Re c 
■ Turbulent, 40 x 80, Re 


' ^ | With cavity (1 : 


0.3 in.) 



Figure 24. Computed laminar and turbulent center- 

line heating distributions at flight conditions. 

Integrated Heat Load 

Carbon-carbon is the forebody heat shield material 
for Genesis. The selection of heat shield material is 
driven by the estimated peak heating rate for the tra- 
jectory. Typically, this peak level occurs at the stagna- 
tion point when the freestream value of pV ' is at its 
maximum for the trajectory. The heating rate at the 
stagnation point is presented in Figure 25 and Table 5 
as a function of time along the nominal trajectory. Also 
presented are laminar and turbulent values at the cavity 
location. Due to the presence of the penetrations, the 
maximum heating rate for Genesis is located down- 
stream of the cavities. The heating rates are two to three 
times as large as the peak stagnation value but still 
within the capabilities of carbon-carbon. 


Figure 23. Comparison of phosphor line cuts with 
turbulent computations (a=0 deg). 


The thickness of the heat shield for a given mate- 
rial is driven by the estimated heat load. The heat load 
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at any point on the body is integrated area under the 
heat pulse as shown in Figure 25. The computed heat 
load at the stagnation point is. The corresponding lami 
nar heat load at the cavity is Q = 9155 J/cnr. In order to 
assess the sensitivity of the heat load at the cavity to the 
prescribed transition point, a number of additional CFD 
solutions were generated (see Figure 25). Assuming the 
How is laminar at the cavities through peak heating 
(t=270 s) and transitions to fully turbulent by peak dy- 
namic pressure (t=280 .v) gives 0 = 16580 J/cnr (es- 
sentially the stagnation value) A heat load based on 
instantaneous transition to turbulence at t = 270 v is 
calculated to be Q = 20,550 J/cnr (25 per cent higher 
than the stagnation value). 

The value of Re ( _, at the cavity is included in Table 
5 and Figure 25. Note that its value at t = 270 s is 194, 
while at t = 280 .s\ Re e = 328. Also note from Figure 21 
that for a 1.0-f'n cavity, the value of w/8 never exceeds 
4 for the nominal trajectory. The wind tunnel data in 
closest proximity to the t = 270 s point on this figure is 
laminar, and there is no fully turbulent data below 
w/8 = 7.0. Given that wind tunnel Hows traditionally 
transition earlier than free flight, a transition to turbu 
lence aft of the cavity at t = 270 .v seems reasonable. 
Further, the assumption of fully turbulent flow' aft of the 
cavities at t = 280 s would appear to be conservative. 
Thus, the current best estimate is that the heating aft of 
the cavities will produce an integrated heat load of the 
same magnitude (Q = 16580 J/cnr) as the stagnation 
value (Q = 16.770 J/cnr). In addition, the inherent con- 
servatism of the wind tunnel, in conjunction with the 
conservative interpretation of the onset of turbulence, 
provides a safety factor for the heating environment 
downstream of the cavities. 



Figure 25. Computed heat load for nominal 
trajectory. 


Table 5. Computed heating parameters along 
nominal trajectory. 


t LvJ 

x/R b = 0. 
q [W/cm 2 \ 


x/R h = 0.7 


q bm [W/cnr] 

1 [W/cnr] 

RCy 

240 

40 

45 



250 

240 

125 



260 

427 

263 

792 

73 

266 

506 

278 

903 

140 

270 

515 

275 

1050 

194 

276 
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219 

968 

274 

280 

335 

159 

783 

328 

284 

232 

118 

498 

349 

290 

111 

55 

155 

359 

300 

25 

10 




CONCLUDING REMARKS 

Cavities in the forebody heat shield of reentry cap- 
sules can produce heat spikes whose magnitudes are 
several times larger than the associated stagnation val- 
ues. Such a cavity can also induce transition to turbu- 
lent flow r downstream. 

For blunt bodies, radial location of cavities has 
minimal impact on downstream heating levels because 
of limited running length. The radial location has a 
larger impact on the strength of the heat pulse, with 
stronger pulses occurring at the foremost positions. 

If a cavity is small enough, no discernable heat 
pulse or turbulent transitioning is produced. Thus, from 
the standpoint of heating environment, multiple smaller 
penetrations are more desirable than fewer large holes. 

The anticipated heating rates for the Genesis SRC 
fall within the capabilities of the carbon-carbon heat 
shield material. Based on the size and location of the 
forebody cavities, the traditional approach of sizing the 
heat shield to the expected stagnation point heat load 
should be adequate for the expected heating environ- 
ment aft of the cavities. 
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